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In large photosynthetic chromophore-protein complexes not all chromophores are coupled strongly, and thus 
the situation is well described by formation of delocalized states in certain domains of strongly coupled chro- 
mophores. In order to describe excitation energy transfer among different domains without performing extensive 
numerical calculations, one of the most popular techniques is a generalization of Forster theory to multichro- 
mophoric aggregates (generalized Forster theory) proposed by Sumi [J. Phys. Chem. B 103, 252 (1999)] and 
Scholes and Fleming [J. Phys. Chem. B 104, 1854 (2000)]. The aim of this paper is twofold. In the first place, 
by means of analytic continuation and a time convolutionless quantum master equation approach, a theory of 
emission lineshape of multichromophoric systems or molecular aggregates is proposed. In the second place, 
a comprehensive framework that allows for a clear, compact and effective study of the multichromophoric ap- 
proach in the full general version proposed by Jang, Newton and Silbey [Phys. Rev. Lett. 92, 218301 (2004)] is 
developed. 

We apply the present theory to simple paradigmatic systems and we show on one hand the effectiveness 
of time-convolutionless techniques in deriving lineshape operators and on the other hand we show how the 
multichromophoric approach can give significant improvements in the determination of energy transfer rates in 
particular when the systems under study are not the purely Forster regime. The presented scheme allows for 
an effective implementation of the multichromophoric Forster approach which may be of use for simulating 
energy transfer dynamics in large photosynthetic aggregates, for which massive computational resources are 
usually required. Furthermore, our method allows for a systematic comparison of multichromophoric Foster 
and generalized Forster theories and for a clear understanding of their respective limits of validity. 



1. INTRODUCTION 



The refinement of laser spectroscopic techniques in past 
few decades has dramatically improved the possibility of 
studying those systems that are at the basis of the presence of 
Life on Earth^ Indeed, the main photosynthetic chromophore- 
protein complexes have been in past decades structurally char- 
acterized at the molecular and atomic level and their abil- 
ity of realizing almost perfect energy transfer processes is 
currently fostering huge research efforts aimed at unveiling 
the fundamental mechanisms at the basis of their remarkably 
high quantum efficiency^ More recently, the presence of 
exquisitely quantum mechanical effects in the energy trans- 
fer process has been revealed in different photosynthetic com- 
plexes even at physiological temperatures^— The observa- 
tions of long-lived electronic quantum coherence in photo- 
synthetic excitation energy transfer (EET) raised questions 
about the role of the protein environment in protecting this 
coherence and the significance of the quantum coherence in 
light harvesting efficiency. In order to elucidate the origin of 
the long-lived quantum coherence and its functional roles in 
EET processes, much theoretical efforts have been devoted 
to construct theories to describe quantum dynamics of pho- 
tosynthetic EET»i^— Generally, quantum dynamic simula- 
tions of large and complex molecular systems require mas- 
sive computer resources, and it is sometimes difficult to ob- 
tain physical insights into them even if such extensive sim- 
ulations are doable. In large photosynthetic chromophore- 
protein complexes, however, not all chromophores are cou- 
pled strongly but the situation is well described by formation 



of delocalized states in certain domains of strongly coupled 
chromophores. The standard Redfield theory— or the mod- 
ified Redfield theory 6 ^— are employed to describe energy 
relaxation within the domains; however, a theory is needed 
for describing excitation energy transfer between the exciton 
states in different domains. For this purpose, Forster theory 63 
has been extended to describe this transfer, namely the gen- 
eralized Forster theory (GFT)^— The original theory must 
be modified in order to include the details of the complexes; 
in particular, when the distance between donor and acceptor 
aggregates is comparable with their physical size one has to 
properly account for the interaction between excitonic states 
belonging to different subunits. Subsequently, the theory has 
been extended 6 ^ in order to account the intra-subunits quan- 
tum coherence; within this framework the distinction in differ- 
ent subunits becomes possible when the intra electronic cou- 
plings among the chromophores belonging to a subunit Ji n t ra 
is much greater than the electronic couplings Jj n ter between 
chromophores belonging to different subunits, and when the 
strength of the system-bath interaction measured by the re- 
organization energy A is such that Jinter "C A < Jintra- 
This approach is termed the multichromophoric Forster the- 
ory (MCFT). It may represent a simplified but effective the- 
ory for the description of energy transfer processes taking 
place in large photosynthetic natural and artificial 6 ^— com- 
plexes where the computational efforts required by the full 
quantum simulations may be prohibitive. Despite these poten- 
tial advantages MCFT has been applied in its full generality 
only in very limited cases and in general its original simpli- 
fied version, the generalized Forster theories (GFT), is usually 
preferred^ Furthermore, a general discussion on the limits of 
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validity of MCFT and GFT is still missing, and this is mainly 
due to the difficulty of expressing the theory within a compact 
framework. 

The main goal of the present paper is to provide such a gen- 
eral and effective framework. The present analysis is based on 
analytic continuation and second order perturbative time con- 
volutionless (TCL2) quantum master equation techniques for 
the evaluation of lineshape operators .r&Z^ which are the main 
tools of the theory. The application of these methods allows 
for a variety of fundamental and interesting results that are or- 
ganized as follows: In Sec. [2] we give a brief review of the 
MCFT and of the main problems that one has to face in or- 
der to evaluate the inter-subunits transfer rates. In Section [3] 
we first review the formalism for the derivation of the exact 
master equations for the emission lineshape operators based 
on analytic continuation and we derive (Sec. 13 Al l a general 
detailed balance condition for the multichromophoric case. In 
Sec. 13 Bl we derive the second order perturbation time convo- 
lutionless (TCL2) general formalism of MCFT and we show 
how its limits of validity can be discussed on the basis of phys- 
ical arguments. In Sec. 13 CI we briefly discuss generalized 
Forster theory and we show how one can derive such sim- 
plified description starting from our general framework. In 
Sec. |4]we apply our theory to simple but relevant multichro- 
mophoric systems: on one hand we discuss the adequacy of 
TCL2 techniques in deriving lineshape operators, and on the 
other hand we study the limits of validity of the full TCL2 
MCFT and compare it with an accurate example of GF the- 
ory. Sectionals devoted to concluding remarks. 



2. MULTICHROMOPHORIC TRANSITION RATES: 
PRELIMINARIES 

Multichromophoric Forster theory is a second order per- 
turbation theory with respect to the electronic coupling be- 
tween different aggregates that can be developed starting with 
the Frenkel exciton Hamiltonian of the overall photosynthetic 
complex. In particular, we consider the following form of 
electronic excitation and bath Hamiltonians: 
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where 
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H NM =J2j™\N,n){M,m\ (2.4) 

nm 

model the electronic Hamiltonian of subunit TV and the inter- 
action Hamiltonian between subunits TV and M, respectively. 
We use capital indices for labeling aggregates and lowercase 
letters for labeling chromophores inside the aggregate: |TV, n) 
is the excited state of chromophore n inside aggregate TV, and 



e„ is its site energy (the Franck-Condon energy). In the en- 
ergy eigenstate (exciton) basis the Hamiltonian for the aggre- 
gate TV reads 
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and 



(N,n\e^) is the change of basis matrix. The 
states |ej^) represent the excitons of subunit TV which are 
delocalized amongst the aggregate's chromophores |TV, n). 
Within each subunit the coupling to the bath degrees of free- 
dom is modeled as 

V£ = £ \N, n) (TV, n\ u» = £ V n N (2.6a) 

n n 

H b N = Y,H^ n (2.6b) 



where £ w H b N = £ ? ?w b (Z) [p(£) 2 + g(£) 2 ] /2 model a 
set of harmonic oscillators, p(£) and q(£) are dimension- 
less coordinates, and is the frequency of the phonon 
mode £. The bath operators are defined by u% = 
— huj b (£)d(t;)n where e?(£)^ is the dimensionless dis- 
placement of the equilibrium configuration of the £th phonon 
mode between the ground state and the excited electronic state 
of the nth site in the Nth aggregate. The reorganization en- 
ergy associated with the nth pigment is characterized with the 
displacements as = ^ ^ b (£)[d(£)n f/ 2 - lt is assumed 
that [flf, H%] = [V£, H^} = for TV + M i.e., the baths 
of different subunits are independent. 

The MC approach consists of an electronic excitation hop- 
ping model among aggregates. As shown in Appendix [A] 
MCFT can be derived from the full quantum evolution of the 
exciton-phonon system by: i) applying second order perturba- 
tion techniques with respect to H , based on the assump- 
tion that J^m ^ Jnm > ^> and ii) assuming that the donor ag- 
gregate D equilibrates on time scales which are much shorter 
than those characterizing the transfer to the acceptor A. The 
first approximation permits to obtain a classical-like mas- 
ter equation for the aggregates populations, with multichro- 
mophoric rates. The equilibration assumption allows to sim- 
plify the problem in that the rates are independent on the ac- 
tual donor's initial excited state; and the non-equilibrium dy- 
namical features of the donor subunits are thereby neglected. 
The resultant multichromophoric (MC) rates read 
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In the above expressions we used the following notation for 
Fourier transformed operators f(u>) = J^° dte lut f(t) and 
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we have defined the MC lineshape operators (LSO) as 



diagonal terms: 



A» a ,(t) = (N,a\Tr 

b 



F»,(t) = (N,d\Tr 





N,a'), 
(2.10) 

\N,d'), (2.11) 
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^A« k ,(tM)(e%\ 

k k^k' 

(2.12) 

where A kk ,{t) = (e^ r |A Ar (t)|e^) and and analogously for 
F N (t). Correspondingly one can write k A< ^ D = k^~ g D + 

^ofHiag anc ' tne off-diagonal contribution is meant to take into 
account the complex intra-subunit dynamics. In the litera- 
ture typically only the diagonal terms are considered, and this 
choice leads to the generalized Ftirster rates 



where the properties A N (t)t = A N \-t) and F N (t)^ = 
F N (—t) are satisfied. Trt O indicates the partial trace of an 
operator O over the bath degrees of freedom. In the com- 
pact expressions, Eqs. (12.71 and ( 12.81 a matrix multiplication 
between the lineshape operators and the matrix J AD with ele- 
ments J^f is understood. The two thermal states are defined 



by p, 



and pb 



where /3 is the inverse 



temperature and the normalization is understood. 

The above expressions are the starting point of the MC 
approach. In order to assess its validity there are two or- 
ders of problem to confront. One set of problems is related 
to the theoretical techniques employed for the evaluation of 
the LSOs. Explicit expressions are known for monomer a 80 ' 81 
but, in general cases, Eqs. ( 12.10b and ( 12.1 It can not be ex- 
actly evaluated because of the interaction, mediated by the 
bath, between the aggregate's excitons. In order to proceed 
further one has to employ perturbative techniques with re- 
spect to the system bath interaction, simulating the dynam- 
ics of the operators A N (t) and F N (t) via non-Markovian 
quantum master equations. One of the standard techniques is 
a second-order perturbative time-convolution (TC2) quantum 
master equation. However, this approach has fatal flaws, even 
for monomers, in producing absorption lineshapes in the non- 
Markovian regimes corresponding to natural situations. It 
is thus obviously outside the bounds of possibility to apply the 
TC2 approach to multichromophoric situations. One is there- 
fore forced to look for a different approach, and this naturally 
leads to the application of the second-order perturbative time- 
convolutionless (TCL2) method, which is able to reproduce 
accurate absorptive lineshape of a monomer. 

The other subtle point that has to be taken into account is 
the very difference between absorption and emission LSOs. 
Owing to the equilibration hypothesis of the donor subunit, 
the emission LSO in Eq. ( 12.1 \\ has an "initial state" p^ b ^ 
1 ® p 1 ^ ■ This introduces further complications as the elec- 
tronic degrees of freedom are correlated with the bath. In 
the TCL2 approach to the LSOs developed in the present pa- 
per this issue is treated in a simple way exploiting an identity 
which relates absorption and emission lineshape: F(t) is ob- 
tained from A(t) owing to the analytical continuation identity 
F(t) oc A(t — ihj3). This point will be extensively discussed 
in Sec. [3] 

The second set of problems that one has to face in order 
to asses the limits of validity and the accuracy of the MC ap- 
proach is the comparison with the alternative, simpler tech- 
niques that are used in the vast majority of the literature. In 
general, the basis used for the expression of k A ^ D is irrele- 
vant, however one can express the multichromophoric LSOs 
in Eqs. (12.10b and ( 12.1 It in the exciton basis and see that 
the distinctive feature of MCFT lies in the presence of off- 
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where J A *' D « = E ad U A k J A d 
ceptor and donor unitary operators that allow to diagonalize 
the respective electronic Hamiltonians. The usual Forster the- 
ory is included as a particular case: in Ref&L it has been shown 
that the multichromophoric rate k A< ^ D reduces to the Forster 
rate k **~ Dl when only two specific optically allowed exci- 
ton states Ak and Dn have non-vanishing average electrostatic 
interactions J Ak > De while all the other j A h-D m ^ (h,m) ^ 
(k, t) are vanishingly small, e.g. in the reaction center of pur- 
ple bacteria. Another example of diagonal theory has been 
employed for describing the rapid excitation energy transfer 
between two circular aggregates, B800 and B850, in Light 
Harvesting Complex 2 (LH2) of purple bacteria^i^i 

While the above approximation is widely used, the frame- 
work developed in our paper will allow us to investigate the 
role of the off-diagonal terms in determining the overall trans- 
fer rate. In particular, our benchmark technique will be a diag- 
onal approach based on TCL equations that is in many aspects 
more accurate than the ones used in the literature (Sec J3 C) . 



3. EVALUATION OF THE LINESHAPE OPERATOR 

In this Section we lay down the general theory for the eval- 
uation of the LSOs. We derive the exact expressions for the 
lineshapes operators and this will be the starting point for the 
subsequent derivation and study of the second order approxi- 
mations in the TCL picture. 

Let us consider the absorption lineshape and, in order to 
simplify the notation, let us remove the index of the aggre- 
gate. A master equation for the lineshape operator has been 
derived in Ref^ 8 - using projection operator perturbative tech- 
niques. Here we derive formally exact expressions exploiting 
the Gaussian nature of the bath fluctuations. Setting x A (t) = 
e -iH cb t/h e iH b t/h tnen m tne " aS y mme tnc interaction 

picture" defined by xt (*) = e i{ - H+H ^ t / h X A {t) e-' lHbt / h the 
following equation of motion is satisfied 



xf{t) = --v eb {t) x A (t) 



yielding 
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where V eb (t) = e i(H+H b )t/H e -i(H+H b )t/H and r+ is the 
standard time ordered product, i.e. the Dyson series. It is 
important to stress that Eq. (I3.lt is not a Liouville equations, 
as the operator V e b(t) acts only to the right. The absorption 
LSO follows from Eq. (13. Il l by tracing out the bath degrees of 
freedom. When the system-bath coupling in Eq. d2.6l l is con- 
sidered, the Gaussian property of the phonon bath yields the 
operators u n to satisfy the Wick theorem. Hence, proceeding 
along the same lines of Ref3 we find 

Al(t) = T+e~T? !° ^ K ^ C aa'(^-s')V a (s) V a ,(«')^( ), 

(3.3) 

where C aa ' if — s) = Tr [pb u a {t) u a ' (s)] is the bath correla- 
tion function. Despite not being essential, in this paper, we 
assume that the bath acts independently on the different chro- 
mophores as usual, i.e. 

C aa >{t) = C{t)5 aa ,. (3.4) 

The bath correlation function can be expressed in terms of the 
spectral density J(u) as 

c «> = ^ <«> 

where J(w) has been extended to negative frequencies, 
J (—(d) = —J(cd), The spectral density J(uj) and the re- 
organization energy A are related with each other via A = 
/„ du> J(tu)/ (ituj). The absorption LSO is then 

A(t) = e-T^'T+e-^^o ds fo ds'C(s-s')V a (s)V a (s')_ 

(3.6) 

From Eq. ( 13.3b we obtain an equation of motion for the LSO 
that, in the interaction picture, reads 

d /"* 

-A^t) = J dsT+ [K(t,s)A l (t)] , (3.7) 
where the dissipation kernel is 

K(t, ») = -ic(i-j)^ V a (t) V a (s) . (3.8) 

a 

It is important to stress that Eq. (13.7) is formally exact pro- 
vided that the expansions in Eqs. ( 13.61 and ( 13.81 ) are inserted 
in the time ordered integral. The latter however mixes the op- 
erators entering in Eqs. (13.6) and (13.8) making Eq. (13.7) not 
suitable for practical calculations. 

Trying to proceed similarly for deriving the emission line- 
shape operator leads to complications as the Wick theo- 
rem cannot be applied straightforwardly due to excitation- 
environment correlations corresponding to the Stokes shift in 
X F (0) = F(0). A possible solution to this problem would 
be to consider the equilibrated state p e b as a state resulting 
from an initially factorized state Pf(U) = p e (U) 53 Pb(U) 
corresponding to the Franck-Condon transition and letting 
ti — > — oo. Such an approach can be analyzed with the help of 
time-evolutions of quantum master equations. However, here 



we consider a much simplified description. Our approach con- 
sists in exploiting the identity 

F(t) =A(t-ih/3)/Z with Z = Tr A(-ihj3), (3.9) 

which can be derived straightforwardly from Eqs. (12.10) and 
(12. lib . The analytic continuation identity permits to tackle 
the unfactorized initial condition in a simple way. In the fre- 
quency domain the above identity reads 

F{lo) = A(u) with Z = dw Tri(w) 

Z J 2n 

(3.10) 

and was used also by Sumi in Ref.— for expressing the emis- 
sion lineshape in terms of the absorption lineshape. 

When the aggregate N consists of a single monomer the 
above equations, Eqs. ( 13.6b and (13. 9b can be evaluated exactly 
without any approximations: 

i m o»W=e"' £i/S e- sW , (3.11) 

FmonW = e -*C«-2Wft e -ff(tr (3.12) 

with 

g(t) = i / da f ds'C(s') (3.13) 
^ Jo Jo 

being the line-broadening function.— The emission lineshape 
directly follows from the absorption lineshape and Eq. (13.9b . 
Indeed, it is straightforward to show 

g(t-ihj3)=g(t)*-^2Xt-j3X. (3.14) 

In the subsequent section we derive approximated expres- 
sion for the lineshapes by making use of lowest meaningful 
order expansion, the second order, of the TCL approach. The 
latter allow to account for all specific cases of system-bath 
interactions and their performances will be compared in the 
relevant situations in Sec. |4] 

A. Multichromophoric detailed balance condition 

Before analyzing the TCL derivation of LSO we discuss 
how, owing to the analytical continuation identity described 
above, the detailed balance condition (DBC) can be general- 
ized to the multichromophoric scenarios. In the standard DBC 
the ratio k Di ~ A /k A< ~ D only depends on the temperature and 
on the energy difference between the donor and the accep- 
tor. On the other hand, it is possible to show that when the 
donor or the acceptor are composed of aggregates, DBC de- 
pends also on the details of aggregate, as the intra-aggregate 
electronic interactions, and on the parameters of the bath. The 
multichromophoric detailed balance condition naturally fol- 
lows by inserting the analytic continuation identity in the fre- 
quency domain, Eq. ( |3.10b , into the definition of the MC rate 
in Eq. ( |2.8b . It is straightforward to show that 

k A ^ D _Z^_ Tr[e-^] 

kD^A ~ z d = Xr[e -/3fff b ] (3 ' 15) 
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where H^ b , for N = A, D, is the Hamiltonian in Eq. ( 12.2b of 
the aggregate N which models both the electronic interaction 
amongst chromophores and the coupling with the environ- 
ment. The latter identity in Eq. ( I3.151 l follows by using the re- 
lation between the partition functions Z N defined in Eq. ( 13.9b 
and (13.10b , together with the definition of A N f-i%P). As for 
the general features, the above equation does not depend on 
the electronic interaction amongst the aggregates: it is given 
by the ratio of the partition function of the two aggregates so 
that it relies only on the details of the equilibrated state of 
the single subunits. The latter however, as discussed in Sec- 
tion [32J can be different from the Gibbs state, especially for 
higher values of A, and depends also in a non-trivial way on 
the parameters of the bath. 



B. Time-convolutionless formalism 

In this section we present the main result of this paper. As 
described below the application of TCL techniques allows for 
a neat formal description of the multichromophoric approach 
that is on one hand compact and general and on the other hand 
it allows to pinpoint the main physically relevant issues at the 
basis of the approach. Both the equilibrated state of the donor 
subunit and the absorption and emission lineshapes can be de- 
termined with equations that are formally consistent. Further- 
more, one can separate the off-diagonal contribution from the 
diagonal one, and its relevance can be studied in terms of the 
electronic structure of a given subunit and of the specific bath 
properties (the details of the derivation can be found in Ap- 
pendix Q where the general theory of TCL approach for de- 
termining the LSOs is sketched - together with an analogous 
derivation with TC2 method). 

Within a second order time-convolutionless (TCL2) frame- 
work the absorption LSO can be written as (see Appendix O 



A 



TCL 



ft) = a ex (t) ■ aft). 



In Eq. ( 13.16b , a ex (t) is defined as 



-iHt/h -g(t) 



(3.16) 



(3.17) 



whereas aft) is the matrix solution of the following differen- 
tial equation: 



aft) = -0(t)a(t), 
where the matrix O is given by 

O kk , ft) = ie^"' J ds Cfs) E kk , fs), 
S fcfc ,(*)=£s£' fe < (e-W-l), 



(3.18) 



(3.19) 



with u k ",k< = u k >< - u)k> and S^ fe/ = Y, a U akU^ kll U ak ' . 
Similarly, owing to the analytic continuation identity in 
Eqs. dig and ( TTUl l we find 



where f ex ft) is given as 

f ex (t) = e~ m e -«(ff-a*)*/i e -»'W (3.21) 
and fft) is the matrix solution of the differential equation, 



fft) = -Oft - W)f(t) 



with the initial condition 



/(0) = a(-ih(3)/TT:A TCL (-ihf3) 



(3.22) 



(3.23) 



The latter is obtained by numerical integration of Eq. (13.18b 
in the imaginary time. 

Equation ( 13.161 ) states that A TCL ft) is given by the prod- 
uct of a matrix a ex ft) diagonal in the exciton basis and the 
matrix aft) that contains the off-diagonal terms. This decom- 
position allows to separate the different contributes of multi- 
chromphoric lineshapes and to discuss their physical interpre- 
tation. The a ex (t) factor describes the lineshapes correspond- 
ing to excitons of the aggregate: the fcth diagonal element is 
the lineshape corresponding to kth exciton with energy hui k . 
Note that each excitonic lineshape can be determined exactly 
without any approximation: the interaction with the bath is 
taken into account by means of the broadening function gft). 
On the other hand, the interaction between the excitons, me- 
diated by the bath degrees of freedom, is taken into account 
by the matrix valued quantity aft) whose physical interpreta- 
tion is more intricate. The matrix operator 0(t) depends, on 
one hand, on the relation between the exciton and site basis 
through the tensor S£ . , ; when the resulting matrix S(t) is di- 
agonal so it is aft), and this in particular happens when site 
and exciton basis coincide. On the other hand Oft) depends 
on the relations between typical frequencies of the system, 
LJk",k', and typical frequencies of the bath that determine the 
behavior of Cfs); in particular if U)k,k' are smaller than all the 
typical frequencies of the bath (high-temperature, Markovian 
case), then e l " fc fc ' t ~ 1: the operator Off) is essentially null 
on the relevant time scales of the bath and hence aft) is close 
to the identity; in this case the behavior of the aggregate is 
well approximated by the lineshapes of independent excitons. 
When typical frequencies of the aggregate are comparable or 
greater than the bath's ones, then the details of the aggregate 
have to be considered and they can give rise to a modification 
of the diagonal elements or to the occurrence of some off- 
diagonal terms. Relevant examples are discussed in the next 
sections. 

As for the emission lineshape, we finally note that since 
the determination of the equilibrated state FfQ) = /oJ CL = 
Tr?, p e b involves the same formal structure previously de- 
scribed one has: 



TCL 



a ex f-ih/3) ■a(-ih0)/TrA l ^ L {-ihP), (3.24) 



F TCL ft) = f ex ft)-fft), 



where a ex (—ih/3) = e~^ H ~^ and af—ih/3) is the solution 
of af—ihr) = —Of—ihr) af—ihr) where the integration in- 
terval is [0,j3 = 1/fcsT]. Therefore, the same above dis- 
cussion holds. Indeed, the structure of the operator Of—ihr) 
is the same as before: the details of the aggregate Hamilto- 
nian are encoded in the same tensor , , , while the "dynami- 
(3.20) ca j" aS p eC ( S are determined by the comparison of the relevant 
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frequencies of the system u)k"k' with the thermal correlation 
frequency lot = 1//3H. In particular, when LOyiy <C ojt 
e.g., high temperature regime, one has that e~ TUJkk ' « 1 and 
the integration interval is small; therefore a(—ih/3) w 1, the 
off-diagonal term in /oJ CL are vanishingly small and the only 
effect is a renormalization of the exciton energies that in gen- 
eral, and depending on A, can lead to an equilibrated state dif- 
ferent from the Gibbs state. In the opposite case, LJk"k' 3> 
i.e., in particular at low temperature, the integration interval 
grows and the off-diagonal terms in pJ CL come into play. 

The above discussion clearly highlights that within the 
TCL2 picture developed in this paper, each steps of the MC 
approach can be easily analyzed by studying a single matrix 
valued operator O that embeds all the structural and physical 
features of the given system under observation. Within this 
picture all further possible and relevant approximations can 
be discussed in a clear way. 



C. Generalized Forster Theory 

Most of the literature about electronic energy transfer be- 
tween molecular aggregates employ generalized Forster the- 
ory (GFT) for estimating the transition rates between multi- 
chromophoric donors and acceptors. ^ 67 ' 77,83 Sometimes the 
internal exciton dynamics within the aggregates is included 
phenomenologically with the use of Redfield/modified Red- 
field theories!^ 3 . However, similar to the assumptions of 
MCFT, usually it is assumed that the donor has reached its 
equilibrium state p e before the transfer takes place. In the lat- 
ter case, the rate between aggregates is given by the weighted 
sum of the Forster rate k Ak ^ Dj from every exciton Dj of the 
donor to every exciton Ak of the acceptor, and the weights are 
given by the the probability of finding the donor aggregate in 
the exciton Dj. The latter are given by the equilibrated state, 
which in the simplest description is modeled by a Gibbs state, 
although this is not accurate because of system-bath interac- 
tion, which causes the Stokes shifti^l Within the formalism 
developed in this paper the equilibrated state is absorbed into 
the definition of the emission lineshape; in this way one can 
show that GFT rates can be obtained as a limiting case of the 
MCFT rates, provided that the lineshapes are diagonal in the 
exciton basis, Eq. ( 12.131 ). 

Equations d3.16t -( [3~23l are the starting points for imple- 
menting various approximations and thus for deriving differ- 
ent GF expressions in a controlled way. A first way for obtain- 
ing a GF expression is to neglect a(t), f(t) and a(—i/3h) terms 
in the evaluation of the absorption, emission line shapes. This 
approximation is a very simple one: the excitons lineshapes 
can be analytically evaluated, and the equilibrated state is the 
aggregate Gibbs state. 

A second GF expression can be derived following the stan- 
dard approach for evaluating the lineshapes in GFT, which 
is based on the off-diagonal elements of a master equation 
which describes the dynamics within each aggregate.- 9 In this 
case TCL2 approximation is usually employed, although with 
some exceptions^ as it reduces to the exact expression for 
monomers. One can show that the GFT lineshapes can also 



be obtained by implementing a secular approximation in the 
general expression, Eq. ( 13.19) : oscillating parts in Eq. ( IC.ll ) 
are selectively removed and the diagonal and off-diagonal el- 
ements are completely decoupled. In particular, in the exciton 
basis 




(3.25) 



Similar expressions hold for A^\, (t), however, since 
A% F k , (0) = 0, one has that A^ F k , (t) = for fc ^ fc'. In 
other words, the secular approximation decouples the diago- 
nal terms from the off-diagonal ones, forcing the LSO to be 
diagonal in the exciton basis. Equation ( 13.251 ) was obtained 
by Renger and Marcus 79 with different methods; furthermore, 
in order to simplify the numerics, they treat in a Markovian 
way the elements with k ^ fc" in ( 13.251 ). Our goal will be 
to compare the accuracy of MCFT and GFT in describing the 
dynamics within the aggregates; to this aim we consider the 
more general expression (I3.25t without any further approxi- 
mation and in the following we will refer to it as GF-TCL2. 
As for the emission lineshape we implement the analytic con- 
tinuation identity, Eq. ( 13.91 in Eq. (13.251 . In this way, no other 
approximations are necessary in the description as the ana- 
lytic continuation identity automatically considers the effect 
of the stokes shift and of the equilibrated initial state. We note 
however that the equilibrated initial states differs from the one 
obtained with the multichromophoric approach, and this will 
be more evident at low temperatures. 

4. ADEQUACY OF MULTICHROMOPHORIC APPROACH 

A. Time-convolutionless lineshapes 

In order to test the adequacy of the MC-TCL2 ap- 
proach we start by considering simple systems. Our first 
goal is to compare the MC-TCL2 rates with those ob- 
tainable with other relevant theories widely used in the 
literature] 19,22 ' 28 i 29 ' 32,38 ' 39 ' 42,46 ' 52 ' 57 One of such approaches is 
the second-order cumulant time-nonlocal (2CTNL) quantum 
dynamics, which takes into account environmental reorgani- 
zation processes or dynamic Stokes shifts^ It implements 
a hierarchical representation^^ 9 - for numerical calculations. 
The second theory is the time-convolution approach for the 
evaluation of LSOs. Although even in the case of single 
monomer's absorption lineshapes TC2 master equations are 
not able to reproduce the real lineshapes and give unphysical 
results, 22,82 the only multi-chromophoric analysis available in 
the literatur e) 68 ' 90,91 use TC2 approximation for calculating the 
LSO. 

We start by analyzing a simple but relevant model i.e., a 
trimer whose Hamiltonian is the following: 

/200 100 2\ 
H = 100 100 2 (4.1) 
\2 2 0/ 
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where the values are given in cm -1 . Since Jintra = 
100 cm -1 ^> Jjntcr = 2 cm" 1 the multichromophoric aggre- 
gate is clearly formed by a dimer subunit and by a monomer 
weakly coupled with the dimer. The exciton energies in the 
dimer are 



dimer 



261.8 cm" 



tlUJ 



dimer 



38.2 cm" 



(4.2) 



The values used for this example are typical values that can be 
found for example in biological systems. The system-bath in- 
teraction is implemented by choosing a Drude-Lorentz spec- 
tral density^ 



J{uj) = 2A 



7W 



Ld 2 + "f 2 



(4.3) 



where 7" 1 corresponds to the timescales of the bath-induced 
fluctuation-dissipation process. The reorganization energy A 
characterizes the strength of its coupling with the excitons. 
Moreover, in this section, we consider the system to be at 
room temperature, T = 300 K. The electronic coupling be- 
tween the dimer and the monomer ( Jmtcr) is chosen to be 
small in order to justify the second order perturbation theory 
which is one of the assumption in MCFT 

In our discussion we consider the energy transfer occurring 
from the dimer subunit to the monomer. Indeed, the mul- 
tichromophoric donor requires the evaluation of the dimers ' 
emission lineshapes that involve more and relevant approx- 
imations. The back-transfer process from the monomer to 
the dimer is not considered as it is related with the forward 
transfer process via the generalized detailed balance condi- 
tion Eq. ( 13.15b , and thus it does not convey further relevant 
information. 

The multichromophoric rate in Eq. (12. 8t calculated with 
time-convolution-less equations are shown in in Fig. Q] for 
7 = 530 cm" 1 3> oj± (Markovian regime) and 7 = 53 cm -1 
(non-Markovian regime). 

As anticipated, in addition to the rates given by MCFT 
we have also performed an exact full quantum calculation by 
means of 2CTNL approach: the exact time evolution of the 
trimer density matrix is obtained and then the rates are esti- 
mated with a least-square algorithm. Furthermore, the rates 
with time-convolution techniques with or without the use of 
analytic continuations (see Appendix (O) are shown for com- 
parison. 

In the Markovian regime (upper panel Fig. [TJ our simula- 
tions show that the TCL2 treatment gives excellent agreement 
with the accurate full quantum calculations for all values of re- 
organization energy A. On the other hand the TC2 approaches 
are accurate only in the weak coupling regime, while in the 
strong coupling regime they completely fail to reproduce the 
exact data. On the other hand, it is important to emphasize that 
MC-TCL2 approach, despite being a second order approxima- 
tion, works excellently even in the strong coupling regime. 

In the lower panel of Fig. [T]are the reported the rates for 
the case 7 = 53 cm" 1 , a value which is comparable with 
the fluctuations timescales obtained in experiments. Here, the 
Markovian approximation is not justified anymore, as fvy is 
of the same order of the dimers energies, Eq. (14. 21 . As in the 




FIG. 1. Multichromophoric rates vs. A (cm" 1 ) from the dimer to 
the monomer for 7 = 530 cm" 1 (upper plot) and 7 = 53 cm -1 
(lower plot). In the above, the MC curves are evaluated with the 
MCFT, where the lineshapes are obtained with different theories. In 
particular, MC TC2 is obtained from Eqs. i flOt and ( 1631 In MC 
TC2+AC, the absorption is still calculated with Eq. dB.4t while the 
emission is obtained from Eq. JB.6t . 



Markovian case the agreement between the full quantum and 
the MC-TCL2 description is extremely good for all values of 
A. 

The behavior of the MC-TC2 rates is more complex. The 
introduction of an analytic continuation approach in the non- 
Markovian do not lead to significant differences. While the 
rates evaluated with this theories coincide and qualitatively 
reproduce the exact one, a close look to the lineshapes in this 
regime shows that TC2 master equations are not able to give 
the correct physical description of the LSO: even for small 
values of A MC-TC2 lineshapes show two peaks. This phe- 
nomenon, which becomes more evident for increasing A, as 
already discussed in Refi^, is completely unphysical. On the 
other hand, this kind of drawback does not affect the line- 
shapes derived with TCL2 approach for all values of param- 
eters analyzed. However, the results of this section show that 
the MCFT approach succeed in determining the correct trans- 
fer rates, and a second order treatment of the lineshapes, as 
given by TCL2, is sufficient. 

We conclude this section by analyzing the behavior of the 
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A(cm- 1 ) 

FIG. 2. Rescaled rates for the Hamiltonian @3) for T = 300 K and 
7 = 53 cm" and different values of Jmtcr 



multichromophoric approach when the distinction between 
subunits is no longer possible. In particular this is hap- 
pens when the condition Jj n tcr < A is violated. In Fig. [2] 
we show the simulations for the normalized rate kn-^A = 
(h/ J inter) 2 ko-t A for the Hamiltonian in Eq. (14.2b and for 
growing values of Ji n t or = J 13 = Jrs- Owing to the rescaling 
the MC-TCL2 rates are represented by a single curve and they 
are compared with the rates obtained with full quantum simu- 
lation approach. While the MC-TCL2 method is always able 
to capture the rates' qualitative behavior, a quantitative agree- 
ment for large Jintor is obtained only for high values of A. 
The result is expected since for growing values of Ji n t e r the 
main hypotheses on which the multichromophoric approach 
is based are no longer verified, and in particular the dynamics 
within the dimer do not lead to an equilibration of the sub- 
unit on timescales shorter than those characterizing the energy 
transfer to the acceptor. 



B. Multichromophoric vs. generalized Forster 

In order to study the role of the off-diagonal terms and thus 
describe the differences between MCFT and GFT we consider 
this simple model of a trimer 

(El J12 Jl3\ 
J12 E 2 J 23 (4.4) 
Jl3 J 23 ) 

where the first two sites form a strongly connected dimer, 
whereas the coupling between the dimer and the trimer is 
small i.e., J12 ^> J13, J-23- The dimer subunit can be diag- 
onalized by an orthogonal matrix U{9), i.e. a rotation in the 
energy space of the pigments in the dimer with some angle 6. 



In the new basis the Hamiltonian takes the form 

e + S J+\ 
e J (4.5) 
J+ J_ / 

As already mentioned the relation between exciton and site 
basis encoded in U(6) determines the relevance of the off- 
diagonal terms in the LSOs. Indeed, in the exciton basis the 
operator in (13.19) that determines the dynamics of a(t) reads 

" ( " ) _ 4 ^ -a sin 4(9 2a* sin 2 26) ' ( ' 
a = l~e- lSt . (4.7) 

In particular the latter is diagonal when 6 = i.e., the ex- 
citons and site basis coincide and each site interacts indepen- 
dently with its own bath, is again diagonal when 6 = 11 /A 
i.e., when the sites energy are degenerate and the excitons are 
given by the symmetric and antisymmetric superposition of 
site basis states. In this case one can write the dimer-bath in- 
teraction term in the exciton basis in the following way: 

Ui = l(6i + 62) + (Tx(h - 62), (4.8) 

1.2 

where a x = |1)(2| + |2)(1|, and consequently one can re- 
define the baths in terms of the center of mass coordinate 
b + = (61 + fe 2 )/v / 2 and the difference between the baths 
coordinates &_ = (u\ — u-i)jv2. Correspondingly one can 
decompose the dimer's Hamiltonian as H = H + + with 
H+ = sl + l®b + + H b +,H- = 5\e + ){e + \+a x ®b^ + H b - 
and [H+,H-] = i.e. the two baths are independent. Due 
to symmetry of this specific interaction, the baths are unable 
to distinguish the excitons. In particular the a x terms swaps 
the two excitons and the due to the bosonic properties of the 
baths the interaction term in Eq. (13.6b (<7 x e -< *' e+ '' e+ ' t o" a! ) is 
proportional to the identity, therefore no off-diagonal terms 
appears in the LSOs. This fact, due to the bosonic properties 
of the bath and Wick's theorem, is not specific to second or- 
der perturbation but it applies at any order of the perturbation 
treatment. 

Aside from the previous two special cases, in order to de- 
scribe the relevance of off-diagonal terms one can focus on 
Eq. d4.6t and fix the conditions for having maximal off diag- 
onal elements. A necessary (although not sufficient) condi- 
tion for maximizing the role of the off-diagonal terms is to set 
6 = 7r/8. See Eq. ( 14.6b . The other free parameters are (see 
Eq. (14. 5b : the off-set with respect to the monomer's site en- 
ergy e, the difference between the dimer's excitons energies 
5, and the electronic couplings J_, J + . We chose to focus on 
the following dimer Hamiltonian: 

(253 63 \ 

^ 63 126/ K ' 

where 6 = 180 cm -1 and e = 100 cm -1 . Indeed, as already 
discussed in section [3131 the role of the off-diagonal elements 
becomes important when the typical frequency of the system 
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FIG. 3. Emission lineshapes for Hamiltonian, Eq. <4.9t ; T = 300K, 
7 = 53 crrr 1 ; A = 10 cm" 1 (a) and A = 40 cm- 1 (b) 



S/H is comparable with the smallest frequency of the bath. 
In the non-Markovian high-temperature (T = 300 K) regime, 
the latter corresponds to 7 = 53cm- 1 /fi,(10ps" 1 ) < 5. 

In terms of emission lineshapes the presence of the off- 
diagonal terms is revealed in the global MC rate when: i) the 
monomer's absorption lineshape sufficiently overlaps with the 
off-diagonal terms; ii) the effect is not screened by the overlap 
of the monomer's absorption lineshape with lineshape corre- 
sponding to the highest dimer's exciton; this latter point can be 
achieved if the two excitons' lineshapes are sufficiently sepa- 
rated in frequency, and this is determined by difference in en- 
ergy between the excitons 5, and if A is not too high. Indeed, 
the emission lineshapes are shown, for A = 10 and 40 cm -1 
in Fig. |3] the absorption lineshape (not shown) is centered in 
zero. When A grows the lineshapes are more broadened and 
the monomer's absorption lineshape strongly overlaps with 
both diagonal terms, thus reducing the relevance of the off- 
diagonal terms. 

In Fig. [5] we compare the rates obtained with MC-TCL2 
approach for the Hamiltonian, Eq. (14.9b with the ones ob- 
tained with full quantum simulations and with the generalized 
Forster theory described in Eq. (I3.25t (GF-TCL2) where the 
lineshapes have been calculated by imposing the decoupling 
of the diagonal terms. Our results show that the multichro- 
mophoric approach properly reproduces the accurate rates In 
general the transfer rates are very small, as we have cho- 
sen small electronic couplings amongst the dimer and the 
monomer ( J13 = 2 cm -1 , J23 = 0) in order to maintain the 
distinction between the two subunits. The figure of merit able 
to quantify the relevance of the off-diagonal terms is the rate's 
relative error Afc = |&gf — &mc)|/&:mC an d it is plot in the 
inset of Fig. 0] We can see that the error between the MC- 
TCL2 and generalized Forster approach is quite relevant up 
to A w 20 cm- 1 ; for A = 10 cm" 1 , Afc is around 27%. The 
behavior shown by Afc is typical for the Hamiltonians with 
9 « 7r/8 in the high temperature regime: the effects of the 
off-diagonal terms is in general quite relevant in the "true" 
MC regime i.e., when Jj n ter <C A < Jmtra and the second 
order perturbation theory at the basis of the approach holds . 
When A ^> Jintra the system enters in a Forster regime; as de- 
scribed before, in this region the presence off-diagonal terms 
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FIG. 4. Rates for the Hamiltonian in Eq. {43 for T = 300 K and 
7 = 53 cm- 1 . In MC-TCL the LSO absorption lineshape is given 
by Eq. J3.16b . in GF-TCL is given by Eq. d3.25b . In every case, the 
emission lineshape is obtained from the respective absorption line- 
shape owing to the analytic continuation identity, Eq. A3. 9b , The inset 
shows the rates' relative errors Afc for different values of A. 
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FIG. 5. Rates for the Hamiltonian in Eq. (49)! for T = 100 K and 
7 = 53 cm- 1 . In MC-TCL the LSO absorption lineshape is given 
by ( f3T6) , in GF-TCL is given by 0,25b , In every case, the emission 
lineshape is obtained from the respective absorption lineshape owing 
to the analytic continuation identity, Eq. Q.9b . The inset shows the 
rates' relative errors Afc for different values of A. 



of the LSOs is screened due to the broadening of the excitons 
lineshapes, therefore the generalized Forster methods are in 
good agreement with the MC-TCL2. 

We finally examine the effect of the temperature. We use 
the following values 7 = 53 cm- 1 and T = 100 K and 
the Hamiltonian, Eq. ( 14.9b . For these values of temperature 
KbT ps 70 cm" 1 , it is comparable with 7 and the typical 
energies of the system. The results of our simulations in 
figure [5] show that again we have two distinct regions. For 
A < 90 cm- 1 the MC-TCL2 approach is in complete agree- 
ment with the accurate rate, while the GF-TCL2 is far less ac- 
curate; indeed, as shown in the inset of Fig. [5] the relative error 
Afc between the two rates is very large up to A = 100, and it 
reaches its maximum Afc = 65% at physiological values of 
the reorganization energy, A = 20 cm -1 . On the contrary, for 
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FIG. 6. Emission lineshapes for the dimer of Hamiltonian in Eq. l |4.9t 

for 7 = 53 cm -1 and T = 100 K, A = 100 cm" 1 (a) A = 
210 cm" 1 (b) 



very large values of A we see that the best approximation is 
given by a generalized Forster approach. This example high- 
lights the limits of validity of the MC-TCL2 approach. For 
values of reorganization energy typical of biological systems 
and when Ji ntcr <C A < Ji n tra the multichromophoric ap- 
proach is extremely accurate in determining the transfer rates. 
For very high values of A and at low temperature the method 
at becomes less accurate; in these region of the parameters 
the equilibrated state of the dimer starts to significantly differ 
from the actual one; in particular its off-diagonal terms are or- 
ders of magnitude smaller than the correct ones. Furthermore, 
and independently from the equilibrated state used, the off- 
diagonal terms in the emission lineshape equations lead in to 
significant errors in the determination of the overall rate. This 
fact is less relevant as for the back transfer rate (not shown), 
where the agreement of MC-TCL2 and the accurate result is 
fairly good. 

For the sake of completeness we include the emission line- 
shapes as given by MC TCL2 for two values of A where the 
MC TCL2 differs from the accurate rate. 



5. CONCLUSIONS 

The study of energy transfer in large light harvesting com- 
plexes can be simplified when these structures can be di- 
vided into weakly coupled subunits of strongly interacting 
chromophores, i.e., when Jintcr <C A <C Jintra, where 
Jntcr(intra) is the inter (intra) subunit electronic coupling 
and A is the reorganization energy measuring the strength of 
the system-bath interaction. In this context, the multichro- 
mophoric approac h 64,68 was put forward to account for the 
intra-subunits bath mediated coherent interaction in the deter- 
mination of the incoherent inter-subunit energy transfer rates. 
In principle, this approach allows to overcome the problem 
posed by the massive computational resources required by the 
sophisticated techniques available for the accurate simulation 
of the quantum dynamics of the whole complexes. However a 
comprehensive discussion of the approach is still missing. 



In this paper we provide a clear, compact and general pic- 
ture of the multichromophoric approach that allows to study 
its general scope and limits of validity. Our discussion is 
based on the evaluation of lineshape operators with a new ap- 
proach based on analytic continuation and second order time 
convolutionless techniques (TCL2). The latter are found to 
be strictly necessary since the alternative time convolution ap- 
proach, despite being widely used in the literature, has severe 
limitations when real biological conditions are considered. 

The picture developed allows one to identify a Forster like 
contribution that can be determined without approximations, 
and that describes the incoherent hopping amongst excitons 
belonging to different aggregates. On the other hand, the bath 
mediated coherent interaction between the excitons of a given 
subunit is taken into account by an independent matrix valued 
quantity whose off-diagonal terms are peculiar of the multi- 
chromophoric theory. 

In this way, one can immediately identify, on the basis of 
the electronic properties of the subunits and of the typical bath 
timescales, if and which multichromophoric (off-diagonal) 
terms are relevant and must be taken into account, and to 
derive simplified versions of the multichromophoric transfer 
rates that encompass the (diagonal) generalized Forster theo- 
ries used in the literature. 

We have applied our general theory to simple prototypi- 
cal examples of photosynthetic systems composed by a dimer 
weakly coupled with a monomer. We have considered both 
Markovian and non-Markovian regimes and we have com- 
pared our results with those obtainable with other accurate 
quantum dynamic simulations and with meaningful instances 
of generalized Forster approximations that can be derived 
within our approach. Our analysis shows that the TCL2 multi- 
chromophoric theory is able to accurately describe the coher- 
ent subunits dynamics and it can outperform the generalized 
Forster techniques especially when Jintcr <A< Jintra- 
The main limitations of the approach occur in the first place 
when the subdivision in subunits is no longer possible i.e., 
Jintcr ~ A and in the less obvious case of low temperatures 
and extremely high values of reorganization energy. 

In conclusion, our analysis provides a comprehensive and 
general framework in which the multichromophoric approach 
to energy transfer processes between complex photosynthetic 
subunits can be discussed. While our numerical analysis has 
focused on simple trimeric structures, the framework devel- 
oped allows for the study of more complex systems and for 
the judicious identification of the aggregates and regimes for 
which a full multichromophoric picture is necessary. 
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Appendix A: Multichromophoric master equation 

General expressions for the multichromophoric rates have 
been obtained in Refs . 64,68 . However, no general formal 
derivation of how MCFT can result from a full quantum ap- 
proach is available in the literature. Indeed, the proof per- 
formed in Refi££ is strictly valid only for a system composed 
of two aggregates, a donor and an acceptor, and does not 
take into account the back-transfer from the acceptor to the 
donor. In this section we complete our theoretical descrip- 
tion of MCFT by deriving the master equation of a multichro- 
mophoric system using a time-convolutionless projection op- 
erator formalism>2i22 

j t Vp eb (t) = K{t)V P eb{t) + l(t)Qp eb (0), (A.l) 

where Q = 1 — V: the dynamics of the multichromophoric 
populations can be obtained from the complete density ma- 
trix p eb (t) by tracing out the irrelevant degrees of freedom by 
means of the projection operator V. We use the notation of 
Section |2] and set 

Vp eq (t) = J2 P n( t )K®Pb, = \N,n)(N,n\, 

N,n 

P»(t) = Tr[(N,n\p eb (t)\N,n)] . (A.2) 

b 

Using a second order approximation with respect to the elec- 
tronic interaction between different aggregate one can show 
that dA.lt reduces to 

M,m 

+ 2(t)Q Peb (0) , (A3) 
with the time-dependent rates 

*£m (*) =\*tf daW™{t,t-a) , (A.4) 
W™(t, s) = Tr (U- t [J] itf W_ S [J] ® Pb ) , 

where U t [] = e~i tC (-), £(■) = T,n\ h *A and J = 
^2 NM H NM . In the following we show the further approxi- 
mations needed to bring Eq. (|A.3t to the form of a Pauli master 
equation for aggregates, i.e. 

jP N {t) = J2 k NM P M (t) - k MN P N (t), (A.5) 

M 

with P N (t) = ^2 n Pn (t). The inhomogeneous term in 
Eq. (IA.3t can be removed formally only if Qp eq (0) = 0, 
namely only if the excitation is localized into a single chro- 
mophore at t = 0. If on the other hand the excitation is 
initially distributed amongst different excitons, even within a 
single aggregate, then an inhomogeneous term is needed. In 
this paper however we do not discuss in detail such an occur- 
rence, as it is expected to affect only the dynamics of the first 
femtoseconds. 



The stationary rates are obtained by taking the limit t —> 
oo, assuming ergodicityiS The multichromophoric rates then 
arise by summing Eq. ( IA.3t for different n and assuming that 
the resulting rates are independent on m. Indeed, 

W™(a)=\im'Ew™(t,t-8) 

t— >oo * ' 

n 

= Tr{H MN U s [H NM ]pf b [m\) (A.6) 
with the equilibrated state 

p£[m] = [n^®/*]. (A.7) 

If the above state is independent on m, i.e. if the state of 
the aggregate M equilibrates to a state which is independent 
on the initial position of the excitation, then the multichro- 
mophoric rate in Eq. (12.7b follows by taking the limit t —> oo 
in Eq. (IA.4t owing to Eq. dA.6t . and assuming that the bath 
acts independently on the different chromophores. 

Appendix B: Time-convolution formalism 

The time-convolution master equation is one of the stan- 
dard methods for describing non-Markovian open quantum 
systems. It has been derived using projection operator tech- 
niques by Nakajima 9 ^ and Zwanzig 94 and has been applied^ 
for the evaluation of the LSOs in Eqs. ( |2.10t and ( 12. lib . Ow- 
ing to the cumulant expansion formalism^^^ it can be shown 
that a second order time-convolution master equation (TC2) 
for the absorption lineshape operator can be derived from (|3.7t 
by choosing the so-called chronological ordering prescription: 

j t A TC (t) = ~H A TC (t) + jf ds K™(t - s) A TC (s), 

(B.l) 

where KZ9(t - s) = e~T: m K(t, s)e^ s , i.e. 

4)W = ~C(t)Y^\a){a\e-i m \a){a\ . (B.2) 

a 

as V a = \a)(a\ where \a) represent the excited state local- 
ized in the chromophore of site a in the acceptor aggregate A. 
Similarly, 

-A TC (t - im =--HA TC (t- ih/3) 
at n 

+ / dzK7 2 ?(t~ih(3 - z)A TC {z) 
Jo 

so that, owing to Eq. d3.9l ), the emission LSO can be evaluated 
by decomposing the above integral along the contour of Fig. [7] 

as 

j t F TC (t) = -^HF^it) + jf* ds K™{t s) F TC ( S ) + /, 

(B.3) 



12 



I I I I I I I I I I I I I I I I I I I 



FIG. 7. Complex contour C for the evaluation of the emission line- 
shape 



where S(ui) = (hut — H)I(2)(ui)/F(0). The above expres- 
sion elucidate the role of the inhomogeneous term; indeed one 
can see that the real part of 7( 2 ) contributes to the shift of the 
the eigenfrequencies of the system (Stokes shifts), while the 
imaginary part contributes to both the width and the envelope 
of the lineshape. 

It is important to stress that Eqs. (IB.4t and (IB. 5b do not 
reduce to Eq. ( I3.12l i when the aggregate consists of a single 
monomer. The following section introduces a master equation 
for LSOs where such a problem does not subsist. 



Wit) = fo drK™(t-ih(P-T)) 



,-tH 



On the 



where 

other hand, the initial state F(0) = A(—iH(3)/Z is evaluated 
using a second order expansion of Eq. (12.10b . Equations (IB. Il l 
and ( IB. 31 ) were introduced in Ref.— for obtaining the multi- 
chromophoric LSO. They are solved with the Fourier-Laplace 
transform as 



J$f(w) = -29f 



w-H/h-iKf°(u) 

F (2) (0)+/ (2) M 
w-H/H-ik%9(u) 



(B.4) 
(B.5) 



where we used the notation /(w) = J °° dt e l "'/(i), and 
F(2) (0) is second order expansion of the equilibrated state. 
The lineshapes obtained with time-convolution master equa- 
tions take the form of a Lorentzian with frequency-dependent 
damping kernel Kj^(u). The above LSOs are not diagonal in 

the exciton basis in general because of KF£{<jS) and J( 2 ) ( w ) : 
in particular the dissipation kernel K?£{uS) is diagonal in the 
site basis. (See Eq. ( IB. 21 ).) As the strength of both the dis- 
sipation kernel and of the inhomogeneous term depends on 
A, and increases for increasing A, the off-diagonal part of the 
lineshape is expected to give a relevant contribution in the in- 
termediate and strong coupling regime. 

Note that the emission lineshape can be obtained also from 
Eq. (|B-41 > once the analytic continuation identity, Eq. ( 13. 10| >, 
is used. We call 



F, 



(2) 



(B.6) 



such a TC2 expression. The latter is more easy to evaluate 
as the second order equilibrated state and the inhomogeneous 
term dot not enter in the expression. As shown in Sec. 14 Al 
Eqs. dB.5l l and ( 13.101 ) with Eq. ( IB. 41 ) give the same results 
when the second order approximation is justified. 

In Eq. ( 13.121 ), the unfactorized initial condition is the origin 
of the stokes shift. On the other hand, in the time-convolution 
master equation it introduces the inhomogeneous term. The 
role of 7(2) (w) in Eq. ( IB. 51 ) can be better understood by using 
the identity 1 + x (1 - + 0{x 2 ), 



-23 



F(0) 



H 



S{u))/h-iK™{u) 



(B.7) 



Appendix C: Time-convolutionless formalism 

Time-convolutionless master equations can be obtained 
from the generalized cumulant expansion by taking the so 
called partial ordering prescription.— ~— In another widely 
used formalism they are obtained from time-convolution mas- 
ter equations formally resuming the perturbative series^ Us- 
ing a second order time-convolutionless approach for the ab- 
sorption lineshape operator, one can show that Eq. (13.7) re- 
duces to 



dt 



-A 



TCL 



(*) 



-H 



dsK™Hs))A^(t), 



(CI) 



where Kf 2 f h (t - s) = e~ im ' n K{t, s)e iHt / n , i.e. 
K^ L (t) = ~ C(t) \a)(a\ e~i m | a)(a | e *H*. (Q2) 



However, in this framework, the implementation the second 
order perturbation approach for the emission lineshape oper- 
ator results in the appearance of an inhomogeneous term that 
describe the existence of unfactorized initial states, Eq. (12.111 ) 

±F TCL (t) = --HF TCL (t) 
dt h 



dsK^(t 



s)F 



TCL 



(*) 



From this expressions it is clear that, while the absorption 
LSO reduces to Eq. ( 13.121 ) in the monomer limit, the same 
does not take place for the emission lineshape. Indeed, the 
above equation can not reproduce the monomer case 



,.Fmon(t) — 



-{e-2\)-g(ty 



F moa (t) (C.3) 



as, in general, I^ h (t) 7^ (i/h)2XF mon (t). In order to over- 
come this problem the emission lineshape is obtained from 
the absorption lineshape via the analytic continuation identity 
CT5J: 



-F 



TCL 



(*) 



-H 



t—ihfS 



dzK^(z)j F TCL (t). 

(C4) 
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In the above formulation the effect of the unfactorized equi- 
librated state is accounted by a different homogeneous terms. 
This is the prerequisite for obtaining a structure that satisfies 
Eq. (IC.3b in the single monomer case: it is indeed straight- 
forward to show that the solution of Eqs. (IC.U and (1C.41 re- 
duces to Eq. ( 13.12b for monomer aggregates, so that they offer 
a more natural expression for the MC-LSOs. They have been 
derived using the analytic continuation identity Eq. d3.9t and 
the second-order time-convolution-less (TCL2) master equa- 
tion. 

The TCL2 lineshape expression can be better understood 
after some transformations. In the exciton basis the TCL2 
dissipation kernel in Eq. iC.2i reads 



% £ S^,C(*)e-^",*' \e k ){e k ,\ 



w 2 



k.k' ,k'' 



(C5) 



where uj k ",k' = Uk>> - Wk>, Z k ' k , = J2 a U a kU 2 k „ U a k- and U 
is the unitary matrix of the basis transformation. Inserting the 
identity one has that 



E-fc" 



k" 



-1) 

(C.6) 



and it is simple to show that Eq. (13.161 l holds, namely that 



A TCL (t) = a ex {t) ■ o(t), where a ex (t) = e -iHt/h- g {t) 
and ait) is the matrix solution of the differential equation, 
Eq. d3lk 



Appendix D: Correlation function 



For the Drude-Lorentz spectral density it is known that 



Cfc e 



k=0 



where 



Co = HX'y ( cot 



4A 



Vk 



1, 

2irk 



(D.l) 

(D.2) 
(D.3) 



The other quantities such as the line-broadening function g(t) 
or the operator elements 0(t) in Eq. (13.19b are then obtained 
from the above expressions. The results are formally summed 
with a computer algebra system. It is found that they can be 
written as a combination of hypergeometric and polygamma 
functions. The latter are then implemented for the numerical 
simulations. 
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